The Casimir effect: medium and geometry 

Valery N. Marachevsky * 
Department of Theoretical Physics 
Saint- Petersburg State University 
198504 St- Petersburg, Russia. 

November 22, 2012 



Abstract 

Theory of the Casimir effect is presented in several examples. 
Casimir - Polder type formulas, Lifshitz theory and theory of the 
Casimir effect for two gratings separated by a vacuum slit are de- 
rived. Equations for the electromagnetic field in the presence of a 
medium and dispersion are discussed. Casimir effect for systems with 
a layer of 2 + 1 fermions is studied. 



1 Introduction 

In the Casimir effect pQ one typically solves electromagnetic or scalar classical 
boundary problems and finds normal modes of the system. The summation 
over eigenfrequencies of normal modes of the system determines its ground 
state energy via a relation E = tkOi/2. In the present paper this relation 
for the ground state energy is applied to dielectrics or metals separated by 
a vacuum slit, which is the typical system in the Casimir effect. Interaction 
part of the Casimir energy of two or more separated bodies in a vacuum is 
always finite, this part determines Casimir forces measured in experiments. 

Different geometries were studied in the theory of the Casimir effect. In 
most geometries eigenfrequencies Ui are not known explicitly, in this case 
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two main approaches exist to derive the Casimir energy. One commonly 
used approach in the Casimir effect is zeta function technique [2J [3]. An 
alternative technique which proved to be efficient in the Casimir effect is 
the scattering approach, it is discussed in Sec. 3 with emphasis on flat and 
periodic geometries. 

Typically Casimir interaction is considered local in space with a given 
frequency dispersion thus allowing nonlocality of photon interaction in time 
direction. Spatial nonlocality of photon interaction arising due to fermions 
and existence of polarization diagrams in quantum electrodynamics make 
the issue technically more complicated. Spatial dispersion is particularly 
important in metals, it yields e.g. a Debye screening. Graphene is another 
example of 2 + 1 system where spatial dispersion can not be neglected. A 
complete theory which would give a description of 3 + 1 systems with spatial 
dispersion in the Casimir effect is still absent. 

A problem of spatial dispersion is closely related to the high temperature 
behavior of the Casimir force. The problem of the high temperature asymp- 
totics of the free energy between two metals attracts attention of theoreticians 
and experimentalists. Asymptotics of permittivity at small frequencies is im- 
portant at large separations between two metals or high temperatures. Two 
models of metal permittivity, the plasma model and the Drude one, yield 
high temperature asymptotics of free energy differing by a factor of 2, see 
Sec.4. 

Instead of using a selected model of permittivity one can alternatively 
evaluate components of the polarization operator in a medium and write 
equations of motions. Components of the polarization operator can be eval- 
uated from the first principles once one knows properties of quasiparticles 
in the medium. This approach allows one to obtain high temperature be- 
havior of the system from the first principles of quantum electrodynamics 
and quantum field theory (see Sees. 4 and 5). Spatial dispersion of the polar- 
ization operator in graphene is important for correct determination of high 
temperature behavior, see Sec. 5 for details. 

A layout of the paper is the following. Sec. 2 is devoted to derivation of 
the Casimir-Polder interaction jl] of an anisotropic atom with a flat surface 
(perfectly conducting and dielectric surfaces are considered) on the basis of 
an approach developed in Ref.[5]. 

In Sec. 3 the scattering theory approach to periodic systems is developed 
after discussion of the Lifshitz formula for two parallel flat surfaces separated 
by a vacuum slit a [5] . Theory of the Casimir effect for periodic systems with 
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frequency dependent permittivity was developed in [?]. Comparison of theory 
and experiments beyond applicability of the proximity force approximation 
(PEA) was performed for grating geometries in [TIE]- 

In Sec. 4 Maxwell equations in a medium with a polarization operator of 
fermions taken into account explicitly are discussed. Free energy between 
two superconductors is analyzed at large separations between superconduc- 
tors. Importance of nonlocality and spatial dispersion in the Casimir effect 
is emphasized. 

In Sec. 5 fermions in a flat 2 + 1 layer are studied following ideas developed 
in Sec. 4. In particular, derivation of reflection coefficients from a layer with 
2 + 1 fermions is given in Sec. 5 in a novel way, diagrams corresponding to 
the Casimir-Polder energy and the Lifshitz free energy for systems with a 
flat 2 + 1 fermion layer are shown. Graphene is a typical example of such a 
system [9], its properties in the Casimir effect have been studied by different 
authors [10]- [16]. Finally the high temperature asymptotics of a graphene - 
parallel ideal metal system is derived. 

Additional references for further reading are given after Conclusions sec- 
tion. 

We use coordinates x 1 1 x 2 ,x i and x,y,z interchangeably. The units h = 
c = ks = 1 are used throughout the paper. 

2 Casimir-Polder energy 

We review general formalism developed in [5] and then derive the Casimir- 
Polder energy of an anisotropic atom above a perfectly conducting plane [I] 
and energy of an anisotropic atom above a flat dielectric. 

The atom is modeled as a localized electric dipole at the point (x 1 , x 2 , x 3 ) = 
(0,0, a), which is described by the current J^x): 

3 

J ( x ) = d^d^ix^ix^ix 3 - a), (1) 

Ji(x) = -di{t)5(x l )5(x 2 )5(x 3 - a), i = 1, 2, 3. (2) 
The condition of current conservation holds: 
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the expectation value of dipole moments is given by the formula [T7] 

/+oo e _i w (t 1 _j 2 ) 
— a jk (u)du, (3) 

where oijk(oj) for uj > coincides with atomic polarizability, a ground state of 
the atom is represented by a bracket ) . We will use the following expression 
for interaction energy E of the atom with a system C: 



E = lim ( In 

Ti->+oo 



Jexp{iS(A) +ij JAd 4 x}DA 
Jexp{iS{A)}DA 



r .(JD C J) 
lim — % - 



n^+oo 2Ti 

(4) 

T\ is a time interval, Dc is a propagator of photons in the presence of the 
system C. 

Consider Feynman gauge of vector potentials. A free photon propagator 
in this gauge has the form 

pD^(u,r)e-^, (5) 

■oo ^ 7I " 

where 

^ ) (w 1 r) = V M 7H (6) 

and = (1, —1, —1, —1) (Heaviside - Lorentz units are used). The propa- 
gator satisfying perfectly conducting boundary conditions A |^=o = A x \ z=0 = 
A y \ z=0 = and ^-\ z=0 = on a plane located at z = has the form: 

D 00 {u,r) = -e iM 7(47rr) + e iMri /(47rr 1 ) 

D xx {uj, r) = e^ r /(4nr) - e i H ri /(47rr 1 ) (7) 

D w (u,r) = e i|aj| 7(47rr) - e^ ri / '(4titi) 

D zz (w,r) = e i|w| 7(47rr) + e i|w|r 7(47rri), 



where r = a/ (x' — x) 2 + (?/ — y) 2 + (z' — z) 2 and 

ri = a/ (V — a;) 2 + (y' — y) 2 + (z' + z) 2 . An image method yields opposite 
signs for Dirichlet boundary conditions and equal signs for Neumann bound- 
ary conditions in fl7J). 

Performing integration over time, using (l3~j), (j4"l), ([71) and subtracting the 
contribution of the free photon propagator (JH]) (to derive interaction energy 
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of the atom with a plane at z = 0) one obtains after Wick rotation the 
Casimir- Polder result jl] for energy of the atom above a perfectly conducting 
plane: 



[ +00 dcuJ^ .. ./ 2n d 2 n 



x'=x,y—y,z'=z=a 



1 



647T 2 a 3 



i=l 

+00 

dwf (a u (iw) + a 22 (^))(4w 2 a 2 + 2ua + l)e~ 2lJja + (8) 
+ a 3 3(«w)2(2ua + l)e^ aX 







To derive the result for the atom above a flat dielectric surface it is con- 
venient to apply the Weyl formula [18J. It is instructive to solve the problem 
in Feynman gauge as well. Components of the vector potential reflect from 
the surface in the following way : 

A p ->■ r TE A p (9) 

A ->• -r TM A (10) 

Ai ->• -r rM 4 (11) 

A z -> r TM ^, (12) 

where p is a transverse (perpendicular) direction and I is a longitudinal di- 
rection to the vector (k x ,k y ), tte{u, k x , k y ) is a reflection coefficient of the 
transverse electric (TE) wave and Ttm^i k xi k y ) is a reflection coefficient of 
the transverse magnetic (TM) wave. 

The free propagator has the same form ([6]). To make needed decomposi- 
tion of waves it is convenient to use the Weyl formula: 

dkrdkn 




i{k x {x-x')+k y {y-y')+ y /L0 2 -kl-kl(z-z')) OfxOKy 



4vrr (2tt) 2 7 J 2^/cu 2 - k* - k* 



valid for z — z' > 0. The atom is located at the point (x, y, z = a), the surface 
of a dielectric is at z = 0. According to the Weyl formula the propagator 
is decomposed into infinite number of plane waves having the wave vector 
(k x ,k y , y/io' 2 — k 2 — ky). A projection of A x perpendicular to the incidence 
plane is A x cos(8) exp(—ik z z') (9 = aicta.n(k y /k x )). After reflection from 
the surface it is transformed into A x r TE cos(6 l ) exp(ik z z'), and the projection 
of the reflected wave from the surface back to the yields the factor 
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A x vte cos 2 (#) exp(ik z z'). The projection of A x to the incidence plane and 
back to the x-axis yields the factor — A x ttm sin 2 (6 l ) exp(ik z z'). Thus for the 
reflected from the surface vector potential A x one gets: 



I r r+oo e i(k x (x-x')+k y (y-y')+^/u> 2 -kZ-k2(z+z')) 

A * = TTTTo / / dk x dk y == x (14) 

(2tt) 2 77_ do 2^/w 2 -k 2 -k 2 

x (tte cos 2 : (0) - r TM sin 2 (fl)) 

Putting a; = a;', y = y', z = z' = a it is convenient to perform integration 
over 9 which yields ir. Performing the Wick rotation and assuming reflection 
coefficients depend on k = ^Jk 2 + k 2 one obtains the contribution of D xx to 
the Casimir-Polder energy: 

^ ^\ a «H, 2 dkk rTE ^ k ]z^ iU > k) e -W^* (15) 



(2vr) 2 i xxv i 4 v / ^+F 

Similar steps allow one to obtain the following resulting expression for the 
Casimir-Polder energy of an anisotropic atom above a dielectric: 



E =-77TT5 / du{a xx {iu) + a yy (iuj))x 
K 1 ^! Jo 

x dkk u2(yTTE(yiU ; ' ^ - r ™^' ~ r ™( iu > k ) k2 c -2a^Ti? 
Jo ' ' Woo 2 + k 2 

^ r duazz(iu) r ^^^p (i6) 



Note that in the case of isotropic polarizability of the atom in the xy 
plane oej(iuj) the expression (Tl6|) can be rewritten in the following form: 

/+oo r+oo , 

duo I dkk[U.tl(iu)Du{iu, k) + Uf p {iu)D pp {iu, k)+ 

+ ILf z {iu>)D zz (iu,kj), (17) 
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where 



D u (iu, k) 



1 -r TM {icj, k)y/co 2 + k 2 e- 2aVZFT: 



(18) 



(2tt) 2 2u 2 



D pp (iu, k) 



1 r TE (iu, fc)e- 2av/ ^ 



(19) 



(2vr) 2 2a/w 2 + k 2 



D zz (iuj,k) 



1 -r TM {iu,k)k 2 e- 2aVZ *+ F2 



(20) 



(2tt) 2 2w 2 V^TP 



and Iljf (iw) = n%,(iuj) = ^(iup 2 , Uf z (ioo) = a zz {iu)u 2 . 

The formula ifTTj) can also be derived by changing the basis locally from 
the beginning using in every point of the momentum space (k x , k y ) in the 
Weyl formula the following identity for a local basis change in two dimen- 
sions: = &xpfixv ^ypfiyv = ^Ipfilv ^p^pv (&xi &y &nd 6;, 6p QIG two sets 

of orthonormal vectors) , so that the new orthonormal frame is written in co- 
ordinates l,p, and respective components of the photon propagator Du,D pp 
and the polarization operator Uf{, U.^ are written in this local frame in mo- 
mentum space. 

3 Free energy 

Consider two bodies separated by a vacuum slit. To obtain the Casimir en- 
ergy and free energy one has to determine eigenfrequencies of normal modes of 
the electromagnetic field. Eigenfrequencies of normal modes can be summed 
up by making use of the argument principle [19J, [20] , which states: 



where Ui are zeroes and are poles of the function f(cu, (3) inside the con- 
tour of integration, degenerate eigenvalues are summed according to their 
multiplicities, also we assume that eigenfrequencies Ui may depend on con- 
tinuous variables /3. The equation for eigenfrequencies Ui of the corresponding 
problem of classical electrodynamics is 0) = 0. For the Casimir energy 
<f>(u) = u/2. 

In the absence of dissipation (when eigenfrequencies of Maxwell equations 
are real) the contour of integration in u plane in ( I2TT) first passes around 
eigenfrequencies Ui and branch cuts on the positive real frequency 




(21) 
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it is explained in detail e.g. in Ref. |21j. To obtain free energy one has to 
substitute <f>(uj) = Tln(2sinh(a;/2T)) into f)2ip . Then free energy can be 
written as 



F = --JdP J du ln(2 sinh(w/2T))im— In f(co, (3) = 

d/3 / ducoth(u/(2T))lm\nf(u,/3). (22) 



'o 

i p p+co 

= — I d,B I 

Using the condition f(—cu*,(3) = f(u,(3)* one expands the contour of inte- 
gration above the real frequency axis: u G (— oo + is . . . + oo + is) and adds 
to this contour a semicircle integral around zero frequency in the opposite 
direction (this semicircle integral is a contribution of the zero frequency Mat- 
subara term). The contour above the real frequency axis is moved around the 
positive imaginary frequency axis, the poles of coth(a;/(2T)) (with residues 
2T) yield Matsubara frequencies u n = 2irnT, and free energy is rewritten in 
the form 

/oo ' 
d(3j2 hyf(icj n ,fi), (23) 

n=0 

prime means n = term is taken with the coefficient 1/2. (Branch cuts on 
the real frequency axis disappear once one puts perfectly conducting plates 
at large separations. Scattering of electromagnetic waves between these per- 
fectly conducting plates yields eigenfrequencies ui located on positions of 
branch cuts. Eigenfrequencies ui effectively transform into eigenfrequencies 
of scattering states which form branch cuts when perfectly conducting plates 
are moved to spatial infinity.) In the presence of dissipation one can also use 
formulas fl22]), ([23]), see a review (20]. 

Consider plane-plane geometry when two dielectric parallel slabs (slab 1: 
z < 0, slab 2: z > a) are separated by a vacuum slit (0 < z < a) following 
Ref. [7]. In this case TE and TM modes are not coupled. The equation for 
TE eigenfrequencies is f(cUi,k x ,k y ) = 0, where 

f(u, (3 = (k x , k y )) = 1 - r^ Edown (k x , k y , uj)r^ ] Eup (k x , k y , u, a). (24) 

Here r!^ Edown (k X} k y , ui) is the reflection coefficient of a downward plane wave 
which reflects on a dielectric surface of a slab 1 at z — 0, while r TEup (k x , k y , u, L) 
is the reflection coefficient of an upward plane wave which reflects on a di- 
electric surface of a slab 2 at z = a. One can deduce from Maxwell equations 
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that r^ Eup (k x ,k y ,u,L) = r^ Edmvn (k x , k y , u) exp(2ik z a) (r^ ] Edown (k x , k y , u) is 
a reflection coefficient of a downward TE plane wave which reflects from a 
dielectric slab 2 temporarily located at the position of the slab 1, i.e. at 
z < 0). 

After substitution of (1241) into the formula (123!) one gets TE part con- 
tribution to free energy of two parallel plates. Contribution of TM part is 
obtained in full analogy. Free energy has the form 

T = r £ jj <m0 H(1 _ e _ 2ov ^ r a, r « )(1 _ e ^ rM)l 

(25) 

w n = 2nnT are Matsubara frequencies, prime means n = term is taken 
with the coefficient 1/2, r^' E = r!^ Edown ,r^ E = r^ Edown , reflection coefficients 
are evaluated at Matsubara frequencies. 

Reflection coefficients (Fresnel coefficients) for transverse magnetic and 
electric flat waves approaching from vacuum the flat surface of a medium 
described by dielectric permittivity e(u) are well known: 

r TM (u, k) = (26) 
e(uj)kl ' + kl 
k (vi) _ k (mi) 

r T E(u, k) = -A—, 7 — rr , (27) 



where 



k { : l) = Ju,* - kl - kl (28) 



k^ = Je(uj)^-kl-kl (29) 



The Lifshitz result for two parallel semi-infinite dielectrics separated by a vac- 
uum slit is obtained from fl25|) when Fresnel coefficients r TM (iu, k), r TE (icu, k) 
are substituted into fl25l) [61 [22]. 

The limit e(iu) — > +00 corresponds to reflection coefficients of a perfectly 
conducting metal plane: r TM = +l,r TE = —1. In this case one gets from 
(121)1) energy of two parallel perfectly conducting plates at zero temperature, 
the result by Casimir pQ: 

E = III dUjdk - dk V ln n _ e -2a v /^ +fc 2 +fe 2 7r^_ 

777 (2tt) 3 1 ; 720a 3 ' 1 ' 
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Another limit is the high temperature asymptotics (AnTa 3> 1) of the 
free energy for two parallel perfectly conducting plates. It can be obtained 
by evaluating contribution of the zero frequency Matsubara term in f l25|) : 

FU«ra»x-TjJ __]n(l- e V •)--g^r- (31) 

We will continue discussion of high temperature behavior of free energy for 
models of metals with frequency dispersion of permittivity and for graphene 
systems in Sections 4 and 5. 

Consider now the system of two periodic gratings with a coinciding period 
d separated by a vacuum slit (FigJTJ [7] . For this system one has to consider 
reflection of downward and upward waves from a unit cell —ir/d < k x < n/d. 
Imagine we remove the upper grating from the system. The reflection matrix 
of the downward wave is defined as Ridown then. 

The solution of Maxwell equations with a given m (m is an integer number 
determining the Brillouin zone of the downward wave with given u, k x , k y ) for 
longitudinal components of the electromagnetic field outside the corrugated 
region (z > h) may be written by making use of the Rayleigh expansion [23] 
for an incident monochromatic wave : 

E y (x, z, m) =i£f ) exp(ia m x - if3 m z) + 

+oo 

^2 R { n ^lexp(ia n x + i(3 n z), (32) 

n=— oo 

B y (x, z, m) =1^ exp(ia m x - i(3 m z) + 

+0O 

22 -Ri„exp(ia n x + i(3 n z), (33) 

n=— oo 

a n =k x + 27m/d } (3 2 n = uj 2 -k 2 y - o? n . (34) 

This solution is valid outside any periodic structure in x direction, in our 
notations it is valid for z > h. All other field components can be expressed 
in terms of longitudinal components E y , B y by standard formulas in waveg- 
uide theory. This can be done since the factor exp(ik y y) is conserved after 
reflection of the electromagnetic wave from a grating. 

To construct the reflection matrix one has to find Rayleigh expansions 
with the condition iffi = 0, Im = 1 and with the condition Im = 1, Im = 
0, m = —oo . . . + oo, m is an integer number. For actual calculations one 
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Figure 1: Two gratings with a coinciding period d. 



puts m = —J ... J, where J is an upper and — J is the lower limit in the 
sums (132|) , ( 133]) . In the case of reflection from a grating the reflection matrix 
is given by 

Rldown(kxi kyi ^) 



Imagine now that we remove the lower grating from the system (see Figj3]). 
We denote the reflection matrix of the upward wave as R 2u p then. Reflection 
matrices Ridown, R2u P depend on wave vectors of incident waves, parameters 
of gratings and mutual location of the gratings. The equation for normal 
modes states: 



where ipi is an eigenvector describing the normal mode with a frequency Wj. 
Instead of equation fl24|) one obtains: 

f 2 (u, (3 = {k x , ky)) = det(J - Ri d own(k x , k y , u)R 2up (k x , k y , u, L, s)). (37) 

For every k x , k y the solution of f 2 (^i) = yields possible eigenfrequencies uji 
of normal mode solutions of Maxwell equations. 




(35) 
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(36) 
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Figure 2: A fictitious grating for which one evaluates R2down- 
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Figure 3: The upper grating in FigJTJfor which one evaluates R,2 U p, normal and 
lateral displacements from the fictitious grating shown in Figf2] are denoted 
by L and s respectively. 
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Suppose that the reflection matrix R2down f° r a reflection from the ficti- 
tious imaginary grating located as in FigfS] is known in coordinates (x,z). 
Performing a change of coordinates z = —z\ + L, x = X\ — s (s < d) in 
f )32|) . (1331 . it is possible to obtain a matrix R 2up for reflection of upward waves 
from a grating with the same profile turned upside-down, displaced from the 
lower grating by Ax = s, Az = L (see FigJS]). It follows that 

R2up{kxi ky, iu, L, s) = 

= Q*(s)K(k x , k y , iu, L)R 2d own(k x , k y , iu)K(k x , k y , iu, L)Q(s), (38) 

where R2down{k x , k y ,iu) is a reflection matrix of downward waves from the 
grating in the system of coordinates (x, z) depicted on FigfSJ Here K(k x , k y , iu, L) 
is a diagonal 2(2 J + 1) matrix of the form: 

K(k x ,ky,iu,L) = (39) 

with matrix elements e -- L A/ aj2+fc y+( fca:H d - ) 2 on a main diagonal of a matrix G\, 
m = —J . . . J. Due to exponential factors in fl39|) all resulting expressions 
are finite. The lateral translation 2(2J + 1) diagonal matrix Q(s) is defined 
as follows: 

QM = ( G 2 G ° 2 ) , (40) 

with matrix elements e 2mms / d [m = —J . . . J) on a main diagonal of the 
matrix G 2 - 

The summation over eigenfrequencies is performed by making use of the 
formula ( 121 p . which yields the Casimir energy of two parallel gratings on a 
unit surface: 

E = - — / du 

(27T) 3 Jo 

x \ndet {^1 - R ldown (k x ,ky,iu)R 2up (k x ,ky,iu,L,ip) J, (41) 

here ip = 2irs/d, s is a lateral displacement of two gratings. This is an exact 
expression valid at zero temperature for two arbitrary parallel gratings with 
coinciding periods d separated by a vacuum slit. 
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Free energy on a unit surface T in the system of two gratings can be 
written as follows: 

T{L, if) = — y~] / dk y dk x 

x In det (i - Ridown(k x , k y , iu n )R 2u p(k x , k y , iu n , L, (pU , (42) 

here u n = 2imT is a Matsubara frequency. The n = term is multiplied by 
1/2. The formula (T4"2"l) is valid for arbitrary profile and arbitrary dielectric 
permittivity of each grating. 



4 Polarization 

Interaction of the electromagnetic potential with a current J given by 
CQ), © yields formalism with frequency dependent dielectric permittivity 
e{uj). This interaction is local and gauge invariant. Polarization bubbles 
make interaction of photons in a medium essentially nonlocal. Nonlocality 
in time is reflected in the frequency dependent permittivity e(u). More 
general nonlocality leads to spatial dispersion, which by now led to intensive 
discussions j2l] - [30] . 

It is instructive to write equations of motion in the presence of a medium 
with a polarization operator taken into account: 

d v F» v + W V A V = -f, (43) 

j M is an external current, IP^ is a polarization operator. 

In momentum space equations of motion can be written as follows: 

n 0n E n 

div£ + — = k) (44) 
-tu 2 E m + (rotrotE) m + H mn E n = iuj m (u, k) (45) 

Here divE = ik n E n , (rotrotE) m = (k 2 S mn - k m k n )E n . 

The components of the 3 + 1 polarization operator can be evaluated us- 
ing techniques of quantum electrodynamics in condensed matter. General 
structure of the polarization operator resulting from gauge invariance in the 
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presence of a medium is the following [3T] (there is no Lorentz invariance in 
the presence of a medium, only rotation symmetry is conserved): 



n „, v _- r) = ^W) (46) 

IT m >, k, T) = U;2 "" n ! 4 ° ( ^ fc ' T) + (5™ - ^f) Il'>, k, T) (47) 
Due to these properties the equation (1441) can be rewritten as 

i^L)[l-^^)=f{u,k). (48) 



In this approach the longitudinal dielectric permittivity is naturally defined 
as Ei(oj, k) = 1 — n 00 (c<;, k)/k 2 , which coincides with random phase approxi- 
mation (RPA) result for the longitudinal dielectric permittivity [32] . 

In equation (j4"5"l) the part H mn E n can be expanded in powers of k 2 . The 
leading term of this expansion is 

n m "£1 |fcV o = -^^i^Xn^ - iV 2 - (49) 



k 



k m (kE) 



E m - — ^ ) ( etP (w, k, T)| |fcVo - 1) U 



2 



where 



e| (fa;,*,T) = l- I ^ M (50) 

e tr (^>) = l-5>^ (51) 

are frequency, momentum and temperature dependent longitudinal and trans- 
verse dielectric permittivities. 

In the limit of zero temperature T — > the longitudinal and the trans- 
verse dielectric permittivities coincide for \k\ — > 0, £i(co, k, T) La^. T _^ = 

e tr (uj, k,T)\,^ 0T ^ Q = s(uj), and the term ( 1491) results in 

n™^| |fcVOiT ^ = -(e( w ) - 1)£"\ (52) 
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The term (1521) is a standard one in the theory of a medium without spatial 
dispersion. 

Neglecting dependence on k in the polarization operator we get the stan- 
dard equation for propagation of waves in a dielectric medium with frequency 
dependent dielectric permittivity e(u): 

- u 2 e{u)E m + (rotrot£) m = iuj m {u, k), (53) 

The equation in a medium has the same form (|53|) in the coordinate space, 
where the operation rot in ([SSD is a standard one in the coordinate space. 

At T = one gets \im^^ U 00 (u, k)/\k\ 2 = oj 2 Juj 2 already for non- 
relativistic electron gas, u 2 [ = 47m e 2 /m, u p i is a plasma frequency, n is 
a charge density, m is an electron mass. So general scheme outlined in this 
part of the paper naturally yields Lindhard RPA result for si(uj, k) [32J and 
the plasma model of dielectric permittivity in the limit \k\ — > of the Lind- 
hard formula: Ei(uj, k)\^ Q = e(u) = 1 - oj 2 Juj 2 . 

Several comments must be maid about the plasma - Drude high tempera- 
ture problem in the Casimir effect. Permittivity in metals at small frequencies 
is a Drude one (£Drude(u) = 1 — uj^Juiuj + ij)) because the current is pro- 
portional to the electric field E, and as a result the pole of the first order 
is required in the permittivity. Optical data for metal permittivities e(u) at 
small frequencies are in a good agreement with a Drude model of the permit- 
tivity [33] • Due to this it is challenging to explain results of measurements 
of the Casimir force between two metals at T = 300K which were found to 
be in agreement with a plasma model of the permittivity [31]. 

The main difference of plasma - Drude results follows from the zero fre- 
quency TE Matsubara term in free energy. Zero frequency Matsubara term 
determines the asymptotics of the free energy at large separations ArcTa 1. 
The Drude model predicts tte^ — 0,k) =0 (lim w ^ £{iu)u 2 = 0) and 
TTM(iu — 0, k) = 1, while the plasma model predicts tte{}u — 0, k) ^ 
(lim^o £{ioj)u 2 = co 2 ^ and ttm^ = 0, k) = 1. At large separations between 
two metals the zero frequency Matsubara term yields leading contribution to 
free energy, and the asymptotical result for a plasma model of permittivity 
is two times larger than the result for a Drude model of permittivity. 

In superconductors the current is proportional to vector potential (or, 
better to say, it is proportional to — d^, where <fi is the phase of the wave 
function of Cooper pairs in the ground state of a superconductor). To study 
the limit of large separations one needs the zero frequency Matsubara term. 
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Suppose II tr (u; = 0, k, T) = m^iT) + ^{T)k 2 + 0{k A ). In this case we get 



r TM (iuj = 0,k) = 1, (54) 

k (v2) _ k (m2) 

r TE (iu = 0,k)= ' * (55) 

where 

k { : 2) = sjk + fcg, (56) 
fcf 2) = A /mg(T)/(l + 7 (T)) + A;2 + A;2. (57) 

Note these expressions are valid for small k only, nevertheless for a nonzero 
mo(T) such behavior of n tr (oj = 0, k,T) for small & leads to the ideal metal 
asymptotics of free energy at separations AirTa 3> 1: J 7 ~ —T((3)/(87ia 2 ). 

As follows from equation (1431) . results ( 1541 and ( 155]) correspond to zero 
frequency reflection coefficients for the electromagnetic wave reflecting from 
a superconductor (if Drude model of the permittivity is valid for the normal 
component of the two-fluid model of a superconductor, see also Refs. [351136] ). 
Indeed, the photon field becomes massive with a mass mo(T) as follows from 
equation (1431) . the condition Ut r (co = 0, k)\^^ = m^{T) is the principal 
condition characterizing superconductivity. Equation (T45l) and the property 
iuoB = TotE yield the equation for a magnetic field with a mass mo(T), 
Meissner effect immediately follows from the equation for a magnetic field. 
Zero frequency reflection coefficients (|54|) . (155]) coincide with zero frequency 
reflection coefficients for the plasma model of permittivity if one puts formally 
m 2 ) (T)/(l + r y(T)) = uj 2 v In fact, if one considers plasma model description of 
metals at large separations (or high temperatures) one should try to explain 
at the same time why superconductivity is not present in the same system. 

In my opinion, one of the principal questions in the Casimir effect is to 
understand the conditions under which it is possible to use the approximation 
( 152]) and the resulting equation ( 153]) in the Casimir effect, i.e. when it is 
possible to neglect spatial dispersion. Spatial dispersion (i.e. dependence of 
n mn on the frequency u and the wave vector k) should be taken into account 
in metals and it is important in metals indeed, it yields Debye screening, 
which immediately follows from ei(gu = 0,k). Of course, spatial dispersion 
makes transition from momentum space to coordinate space complicated 
from the mathematical point of view. However, it is worth studying the 
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issue since the solution of the problem of high temperature behavior of the 
Casimir force between two metals can hardly be obtained without a detailed 
study of spatial dispersion and properties of the polarization operator at 
finite temperature. 



5 Fermions in a layer 

In this section we find reflection coefficients for transverse electric (TE) and 
transverse magnetic (TM) modes reflecting from the 2 + 1 layer with fermions 
at z — 0. High temperature results for a graphene - ideal metal system are 
discussed. 
Equations 

lead to conditions 



d^ v + 5(z)U vp A p = (58) 



d z A m \ z=+ o — d z A m \ z= ^Q — U mn A n \ z= Q. (59) 

Let's consider the condition 

d A° + diA 1 + d p A p = 0, (60) 

here a direction of the wave vector in the xy plane is along the coordinate 
I, p is a transverse direction. In fact, the condition ( l60l) is quite convenient 
for a description of transverse electric and transverse magnetic modes of the 
propagating electromagnetic wave. 

For a nonzero A p , the condition d p A p = and conditions A\ = A z = Aq = 
describe propagation of the TE electromagnetic wave (the electric field is 
parallel to the surface z = 0) since E p ~ A p . 

For the TE wave we have: 

A p = e iM e ikzZ + r TE e ikl e- ikzZ for z<0 (61) 

A p = e ik * z e m t TE for z >0 (62) 

and 

U pn A n = -U pp (u,k)A p . (63) 

Here k 2 z = u 2 — k 2 , k 2 = k 2 + k 2 . From continuity of potentials at z = one 
gets 1 + r TE = t TE . Now one substitutes (16TT) and (|62|) into (159!) and uses 
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(163)) to obtain: 

n pp (iu;,A;) n pp (2w, k) ( U pp (iu,k) 

r TE Uu, k) = — . = — , 1 , 

2Vco 2 + k 2 - U pp (ico, k) 2Vco 2 + k 2 \ 2^/u 2 + k 2 

(64) 

Conditions A p = A z = 0, koAo = kA\ describe the transverse magnetic 
(TM) wave. This choice of vector potentials describes TM wave since E z ~ 
8 z Aq or B p ~ d z A\. For Aq we have: 

A = e ikl e ikzZ + r Ao e ikl e- ik * z for z < (65) 
A = e ik * z e ikl t Ao for z > 0. (66) 

Two conditions follow from gauge invariance: 

ujU 00 (uj, fc) - kll l0 (uj, k) = 0, 
ull i(u, k) - kUu{u), k) = 0, 

which yield 

u 2 Ii m {uj } k) = k 2 Ii u (uj } k). (67) 

Thus one gets 

2 h,2 r,2 

U 00 A° + U ol A l = -^—U U A = -^IL U A , (68) 



L0 Z UJ Z 



and from (|59|) the condition 



k 2 

2ik z {t Ao + r Ao - 1) = |n„A (l + r Ao ) (69) 

follows. After use of r™(jw, fc) = — r Ao (iu, k) and Wick rotation one obtains 

^/iJ2 _|_ fc2 ( ^/Z^rk 2 \ 1 

r TM (iuj, k) = — — Uu(iuj, k)\l — % — U u (iuj, k) 1 (70) 

It is interesting to see how these results can be interpereted in terms of 
diagrams. Consider first 2 + 1 fermions interacting with an atom. 

In the gauge Aq = the longitudinal part of the free photon propagator 
has the form (i = 1,2): 



nL( . , , hk 3 y/^TFe-Nv^ 
D^iu, k, z) = — — , (71) 
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Figure 4: Typical term in the Casimir-Polder energy between an atom and 
fermions in a 2 + 1 layer. Thick wavy line represents D(a), thin wavy line 
represents D(0). 

the transverse part of the free photon propagator has the form: 



Consider TE part of the formula ( fl71) . Substituting (164j) into (fl71) and using 
(172]) we get: 



The expression f!73|) has a clear diagrammatic representation (one of the 
terms is shown on FigH]). Sum of the terms in round parenthesis in (1751) is 
the sum of RPA diagrams in the two-dimensional layer (note that photon 
propagators do not depend on z there). Taking the trace of (!73|) one gets 
TE contribution to the Casimir-Polder energy of interaction of the atom and 
2 + 1 fermion layer. 

Consider the formula (1251) for two parallel layers with fermions separated 
by a distance a. One gets diagrammatic representation of free energy ex- 
panding the logarithm. The factor 1/n is a standard factor arising in repre- 
sentation of thermodynamic potential in terms of closed loop diagrams with 
n equivalent clusters connected by n photon lines. The cluster is shown on 
Figj5j typical diagrams contributing to the Lifshitz free energy are shown on 
Fig-El Recent discussions of related multiple scattering techniques can be 
found in Refs.jlE], g7J. 




(72) 



np*-Dpp — lip* Dp p (2a) ■ r TE — 

= n£D£(a) (n pp + U pp Dl(0)U pp + . . .) D T pp (a) (73) 
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Figure 5: Cluster in the Lifshitz formula. Thick ellipse represents sum of 
RPA diagrams in a fermion layer. Thick wavy line represents D(a). 




Figure 6: Diagram representation of the Lifshitz free energy for two parallel 
fermion layers 
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Graphene is a typical 2 + 1 system that can be studied along the lines 
described in this part. Quasiparticles in graphene [5] obey a linear dispersion 
law lo = Vpk (vp ~ c/300 is a Fermi velocity, c is a speed of light) at energies 
less than 2 eV. There are N = 4 species of fermions in graphene. The 
polarization operator for 2 + 1 fermions at finite temperature was found in 

unj. 

Define trll = II™. Due to the property (I6T|) and 



. U 2 ^_k 2 
k 2 



trll(zc<;, k) = n o(ic+ k) — — H pp (iu, k) (74) 



one can rewrite reflection coefficients f)64p and f!70p in the form given in jTHj: 

/ t \ ^n 00 . M fc^IIoo + fc 2 trll 

k) = kJIn + W TTE ^ k) = - k 2 U 00 + k 2 (trU-2ik z y (?5) 

For k — > one gets for the zero mass gap and zero chemical potential (a 
is the coupling constant): 

n/ A/"?) 2 /V 2 

trn(iw = 0, k) - U 00 (zu = 0, k) = ^ +.... 

For the ideal metal ttm = 1,?"t.e = — 1 for all frequencies and wave vectors. 
Zero Matsubara TM and TE terms yield following high-temperature behavior 
of free energy (|25p in the graphene - ideal metal system: 

lbira z 
aNvp 

• FDTE = _ 192rf + -- (77) 

Here 

TC(3) - T I - 1 T I fTtt 

— — n = •/ Drude T-s-oo — — id T-s-oo- {<o) 

lbira z 2 
is the high-temperature asymptotics of the metal - metal system with a 
Drude model of permittivity used [37] - |45j, which is equal to one half of 
the high-temperature asymptotics in the metal - metal system with ideal 
boundary conditions (|3T|) or the plasma model of permittivity used |34j . The 
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zero frequency TE Matsubara term is suppressed by a factor aNv 2 F and 
additional power of l/(Ta). 

The high-temperature asymptotics of free energy in a graphene-metal 
system coincides with the Drude high-temperature asymptotics of the metal- 
metal system. Detailed analysis [10] shows that the high temperature be- 
havior in the graphene-metal system takes place already at separations of 
the order of 100 nm at temperature T = 300K, which makes systems with 
graphene very promising for experimental studies of the finite temperature 
Casimir effect. 

The example of this section shows that nonlocality arising due to in- 
teraction of photons via a polarization operator plays a crucial role in the 
high temperature asymptotics of the system. Vacuum effects arising due to 
fermions are essentially nonlocal and lead to spatial dispersion. One of the 
most challenging problems in field theory is to unify nonlocality in momen- 
tum space with boundary problems in coordinate space. Casimir effect is a 
natural area for further developments in this direction. 

6 Conclusions 

Casimir effect is an area of research where one applies and develops methods 
of quantum field theory in the presence of a medium. Throughout the paper 
I discuss several techniques which proved to be efficient in the theory of the 
Casimir effect. 

Casimir-Polder effect is studied in Feynman gauge of vector potentials in 
Sec. 2, the Casimir-Polder energy of an anisotropic atom above a flat surface 
with boundary conditions of the ideal conductor and above a dielectric surface 
is found. Scattering approach is an efficient method in the theory of the 
Casimir effect, it is introduced in Sec. 3 and applied to flat and periodic 
geometries separated by a vacuum slit. 

Properties of the medium determine important characteristics of Casimir 
systems such as the high temperature asymptotics of free energy once one 
knows dielectric permittivity of the medium. An alternative approach is 
to consider properties of the medium through evaluation of components of 
the polarization operator, in the current paper this idea is being developed 
in Sections 4 and 5. In Sec. 5 exceptional finite temperature properties of 
systems with graphene are discussed. 
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Further reading 

A discussion of interaction in the Casimir effect is given in the paper [48J. 
In the review [20] applications of the argument principle to flat geometries 
are considered, dissipation and dispersion are discussed. Papers [21], 09] yield 
detailed derivation of the Lifshitz result for the plasma model of permittivity. 
Finite temperature Casimir effect is reviewed in [50J, see also a review [51J. 
Examples of scattering theory applied to spherical and cylindrical geometries 
can be found in papers [52] - [55] and in Chapter 10 of the book [56]. Casimir 
repulsion due to Chern-Simons term was studied in Refs.[57J [58] , scattering 
of electromagnetic waves on a plane with Chern-Simons term was considered 
in Ref. |59j. Casimir- Polder interaction and macroscopic QED are reviewed 
in [BO! ED E2], scattering approach to Casimir- Polder effect is discussed in 
Ref.[63j. Recent advances in microstructured geometries and Casimir effect 
are summarized in a review [64] . Techniques for evaluation of functional 
determinants are discussed in [65] , methods of heat kernel and zeta function 
are presented in reviews [661 EZ] and books [2], [3], [68], [69] in detail. One 
can find examples of zeta function technique in the Casimir effect problems 
in papers [70] - [75] . 
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